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Preliminary Remarks. 

1. JlHE investigations presented in this paper consist of two parts ; the first offers 
a solution, in a certain qualified sense, of the general linear equation in finite differ- 
ences; and the second will be found to give an almost complete analysis of the 
resolution in series of the general linear differential equation with rational factors. 
The second part is deduced directly from the results of the first, although the 
subjects of which they respectively treat appear to be wholly independent of each 
other. 

With the exception of a few cases capable of solution by partial and artificial me- 
thods, there does not at present exist any mode of solving linear equations in finite 
differences of an order higher than the first ; and with reference to such equations of 
the first order, we are obliged to be content with those insufficient forms of functions 
which are intelligible only when the independent variable is an integer, and which 
may be obtained directly from the equation itself by merely giving to the indepen- 
dent variable its successive integer values. It is in this insufficient and qualified 
sense that the solutions here given are to be taken ; and the first part of the follow- 
ing investigations may be considered as an extension of this form of solution from the 
general equation of the first order to the general equation of the nth order. 

Linear Equations in Finite Differences. 
2. A complete analytical theory of the general equation of the nth order, 

would involve its resolution into a series of equations of the first order of the form 

U A , r W U X _ i = \j w 



x x) 
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from which the complete solution would result, so far as known methods extend, in 
the form 

7/ _ p' p' p' 5 ( ®x \ j_ 4. p(»)pc») p(»)2 ( G r \ 

^~ ri 2# ^Vpp..pUi/ 2 " Vp (w) pK7p^/' 



or in an expanded form, 

^x^^^x "T^x^x-l \* x-\*? x^*x-Z "T'*x-2Xx-"l*x* J x~3 I ••• I C" m* *n+l- • • * * 

J- IT J_P"fi" -l-P" P'(T 4-P" P" P"(T 4- -l-/"P'P" P" 

I • • • • • 

l~" • • • • • 

JL pW_LP W r W J_P (w) p( w )f^( w ) 1 pW T)(ri) p(»)p(») 1 1 0)p(rc)p(w) p(w) 

T^ff I *# U *-lT r *-i r « U ii;'-2T r *-2 r «-l r « VJ ^-3T"«T t r fl» r w+l'"r« • 

The last column contains, as is well known, the general solution of the original 
equation deprived of the term G, ; and the remaining columns contain the particular 
solution of the original equation. 

3. The investigation the results of which are given in this part, although it actu- 
ally succeeds in solving the original equation, will be found to contribute little or 
nothing to the analytical theory as above explained ; and this arises from the circum- 
stance, that the particular solution, instead of being produced in the separated form 
above written, is produced in an aggregated form, 

G^+ AG*_ 1 +BG,_ a + . . . ; 

and the complementary part of the general solution, instead of being produced as 
above, appears in the form of a sum of the complementary functions above written, 
the constants c', c", &c. being the same in all. But as we shall thus obtain a solution 
with an arbitrary constant, capable of solving the original equation deprived of G w 
we have it in our power to solve the equation completely by reducing the order of 
the equation in successive steps. 

If we take the equation of the first order, 

and arrive at its solution, not properly by any analytical method, but by giving to x 
the successive values 

x 9 #— 1, X-—2 ra-f-1, m, m—- 1, 

and eliminating all values of u x between the first and the last, we have 



• . • • 



V* ** #-!• • ••* mi^m-ir V* v^m+lj^ml (,*#•• " m-\-2)^ J m+l"T' ••• T *x*x~l^x~2 m T'±xl*x-.i "T"G^ 



FINITE DIFFERENCES AND LINEAR DIFFERENTIAL EQUATIONS. 263 



— (P P \Lt 4- — 4-- m+1 4- A — - 



#-2 . VT^_i f VJ X 






p n p p \ v> p 



which is the solution ordinarily given, though it is arrived at by a process somewhat 
less coarse than the above. 

By applying the same process of successive elimination to the general linear equa- 
tion (1.), and carefully observing the law by which the entrance of the factors 
P^ Q^ R,, &c. is governed, it will be found that an exactly similar solution may be 
found in the form 



U x — S 1 **•••* w( ^( p p ) l C ) i) 



or 

^^ = ^ T ^^ j a«-- i T ^ X^tt^x-ij^x—z i £ \** ,M "#-2/ , "#-3T ,, T^ V x #»» •*• w+i/ •*■*»» l ^ V x # ,,,x »i/ 

where ^ denotes the sum of a set of distributive operations V 15 V 2 , V 3 , . . . V w _! not of 
a strictly algebraical character, which are capable of being performed only upon fac- 
torial expressions containing consecutive values of P x , and which have the following 
significations. V! denotes one operation of this character, signifying that the factor 
Pm-i Pm is changed into Q m as often as it occurs, any term in which it does not occur 
disappearing, and the sum of the terms thus obtained being the result of the opera- 
tion ; so that, for example, 

V 1 { r x ^ 2 ^x- i x #/ == *°ix-- 1 x x \ * a— 2^1x9 V i \r#^. 2 a #„ i Y x ) — U, 

* l\Xx~3*x-2±x-lL x ) == W^_2* a?-l x a? I "*-3y»-l"* I ^-3^-2^t^5 ▼ 1 V * *-3 *x-2*x~- l*x) =Z ^Wa-2^ix9 

&c. &c. 

Again, V 2 denotes another operation of a similar character, signifying that the 
factor P m _ 2 Pm-iPm is changed intoR w as often as it occurs, the result of the operation 
being as before the sum of the terms ; so that, for example, we have 

V (P P P^=R 

v 2V X X-2 L x-l*- x) XX X) 

V /p P P P P^- P P p -LP P P ip p p 

v 2\ x x— 4 X x— 3 X #--2 x x—l x x) "— "* iV *'~2 x #-~l x x^ x #~~4 xl/ #— I 1 a? |^ A ^? — 4 X #— 3 XV £'J 

* 2\ x #—5 * ' ,X W *\x— 3 X #— 2* #— l x #""+" X #— 5""**— 2* «— l x # I *X—$*x— 4 X \«'-1 X <3?"T" X #— 5 X «t?-4 x A'— 3*V»J 

^2\^— 5* • ,x ^?j = ~^ x *-^— s*^^ <^C. etc.; 

V 3 denotes the change in a similar manner of P m _ 3 P m _ 2 P w _ 1 P m into S w ; 



V w __ 2 denotes the change in a similar manner of P m _ w+2 ...P W into W 
and 

V w _! denotes the change in a similar manner of P m _ w + 1 ...P m into Z m . 



m ? 
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It might at first sight be supposed that this process, if successful at all, would give 
the complete complementary solution in the form 

U X ~=^C ] Z \ < ir l . . r x ) -f- C 2 S { r 2 . . r x ) -j- . . ; 

but it will be found that in reality this introduces only one arbitrary constant. 

It now remains to place these expressions under a properly algebraic form, and to 
verify the result ; and in order to do this we must express the general term e v (P x ...]? p+1 ) 
in terms of the factors of the original equation ; and afterwards give to p, which in 
the first instance is regarded as a constant, the successive values required for form- 
ing the several terms of the solution. 

This may be done as follows :— 

tt > ^°*w ^x ^x V* x ^^ i £*x ^^ 

LiGt p p- — q X ) p p p — t^ p p p p~ — s a9 . . . p p- — w a9 ana p p- — 2^. 

x x—l*-x x x— 2 X x— l x x x x— 3 X x— 2 L x— l x x x x— n+2»* • x x x x— n+l • • • x x 

Then it is easily seen that 

V i\ x x' • •*p + l) :== *x" • ^p + l\(Jp+2l~ tfp+3* • • • I $xJ * x* • • *p + l^p + l$x+l) 

* 2\*x" m *p + l) == x x . . . ±p + i\'Pp+3~T"Vp+4 : "!"• '*\'x 5 * #••• *p+l^ f p + 2' x+l) 

and generally 

V\r x . . . r p+l ) = r x . . . r p+1 { 2^+!^+! + ^>+2^+i H~ Z p+3 S x+ i-f- . . . + 2, p+n _ 2 W x+l -J- 2 p+n _ 1 Z x+l } . 

To find V' (P X ..P P+1 ), it will be convenient to proceed by steps, beginning with a 

small number of terms. 
Thus 

»1V* #• '**x-3) = — \ x ^••• x x-3j\^xlQx-l I ^x-2Jy 2 * i\*# % " x~3/ x " * x-3\9.xQx-2J 

V 1 (P^...P^_ 4 ) = (P^..P^_4)(^+^_ 1 + ^- 2 + ^-3)^ ^ V l (^- * - P ^- 4 ) = P^ • • .^^-4(^(^-2 + ^-3) + ^-l^-s) 

* l\*X' •* "*-5/ == V x A'' • •*x-5j {tfxltfx-l I 5 r < ^-2'T"5 r A?-3'T"5 r ^-4J5 

T| 

and generally, V 1 (P,...P P+1 ) being P,.. JP p+1 (2* +l g, +1 ), we have 



Similarly, it will be seen that 

2 3 * A* #****p+i/ ==x #' •• x ^+i^+5(^+i^jo+3C^+i^jo+i5'a'+i < /)5 ana so on< 
In like manner we shall find 

~v 2 (P p \~~p p y* ( r y?~ 2 r \ 
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2 3 * 2\ #' * p+i)"" #'* ' p+ i "p+8y# +i p+$v #+i^ , j9+2^Wi/j? &flQ so on, 

v x being that term of the series q x r x s x9 &c. which corresponds to the operation V m , as 
q x corresponds to V 1? r x to V a , &c. 

^V m (P^..P p+1 )=P^...P p+1 2^ +m ^ +1 (4 +1 2^+^+i) ? (4 corresponding to V,) 5 

W * m\*x* • • ^p+l)~~ ^x* • • ^p + l^p+m+2l+2\^x+l^p+m+l+l\^x+l^p + m^x + l))} &C. &C. 

Finally, if N^+N^+NJ jjP +- represent a series in which 
and generally 

^ffjjp ^+2a-lV^+l 1 ^ l ^-l, j p/ I ^p + 3a-l\' x+l L ^ x-2, p} "T" ^p + ia-l^x + l 1 ^ x-3, p) HT " " 

\ £p+{n-l)a-\\W x +\™x-n+2, p) I ^p+na-l\^x+l^ x-n+1, p) 9 

then, if this series be called N^, we shall find that b%P x ...V p+1 ) = (P x ...P p+1 )T^ X)P . 

It may be here observed that the number of the terms of the series Nl )i ,+N^+.. 

l 1 

cannot exceed ^{x—p)-\-l when w—p is even, and ^(^"-"P+l) when x—p is odd ; and 

since p has the successive values x—2, x—3, &c., it is always known whether x~ p be 
odd or even. The number of terms may be less ; for if Q, be zero, the number of 

terms would be the next whole number above i(^— p), &c. 
That the equation 

W / X := -C\*W** * p + l)™x,p 

is a complementary solution of the original equation, or in other words, a solution of 
the original equation wanting the term G„ may be directly verified ; for we have 

* — C\*x* • * p+\) ^*^* x-\, p' 

Now AK. 1|P =0, 

&N x __ liP =q x -\-v x -\-s x -\- ,.-\-z x ^ 

AN'" =o N" 4-r N" -I- v N" 4- 4-?/) N" 4-^ N" 

*-*-^o?— 1, p ^ A ^-2^T' x x ^x-3, p\°x L ^x-4,pl • • • T^r^^-w+l^TV^jy-M,^ 

and generally 

AN (a+1) — n \r (a) -4-r N (fl ° J- c N (c ° -L J_™ lV (a) J_ cy tm (0) 

A - vl ^^-l,p !/^ 1 ^ l ^-2,j9 I ' x 1 ^* x-3, pHT *x^ x-4, pf ••• "r^^^-n+l,^ TV"«-«, 

whence 

^^x-l,p = qx™x-2,p^Fx™x-3,p'~\~SxNx-4 : ,p"T' •• "r^^'*-n+l,p"r^N^_ Wj 

and 

^ ±x U x- 1 === y A?-2~T* ^x U X-3 \ ^ A?-4 ~r • • "T" W A?-n+ j + /j x u x _ n , 
MDCCCL. 2 M 



J?5 



p? 
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In this part of the solution, I apprehend that it will not generally be necessary to 
have regard to the lower limits, since p may have any constant value, and that value 
may be taken which is most convenient in each case; and if we make p=—ac, all 
the lower limits will have this value. 

This part of the solution would then be 

M*=cP,....P 1 {l+A,+B,+C„+..}, 
where 

A#= Ztyx+i + -^Vf i T ^x+i "T* • • "T* 2fWx+ii ZZx+i \ 

and generally each term of the part within the brackets is formed from the preceding 
term thus : first change x into x—l, multiply by q x+u and effect a summation ; then 
change (in such preceding term) a: into x—2, multiply by r x+19 and effect a summation; 
and so on until lastly we change x into x—n-\-l, multiply by s, +1 , and effect a sum- 
mation ; and the sum of the parts thus obtained is the next term. 

The verification of the particular solution is easily derived from the above ; but it 
rests on the assumption that the algebraic value above given for e v (P x ..P p+1 ) is correct; 
to which therefore particular attention is directed. 

The equation 

considered as a solution of the original equation wanting the term G x9 implies that 
Now in order to verify the particular solution, 

U x ==L ^ y x\ *a^*a-i I S \*x*x-l)* Jf x-2~Tt ("** *-i-^*-2y^ jr jif-3 l • • + 2 \*x'*± x _ p + l )KJ x _ p -f~...) 

it is only requisite that 

£ v (r x ..r a ,_ p+l ) = r a e (P#_i..P#-p+i ) + **#£ (r #_ 2 ..x x-p+i)\"^x^ [r x -. 3 .*l x ^. p+ J + .. + Zi^g (P^-»».P^. / , +1 ) ; 

and this is true, for it is identical in form with the equation last above given, not- 
withstanding the occurrence of x in the lowest value of P, for this lowest value re- 
mains the same throughout the expression. 

7. If P,=0, the expressions e v (P x ...l? w „ m ) reduce themselves to those terms which 
do not contain any value of P. It would not be difficult to determine generally 
what terms these are; but probably the most convenient general method of arriving 
at the solution in an algebraical form would be to make P, equal to a constant b, and 
to make b equal to zero in the result finally obtained. 

8. The above particular solution of the original equation is in such a form that the 
general term of the indefinite series representing the solution is given in explicit 
terms ; but that general term may be represented in an implicit form, which perhaps 
is more convenient for practical use. 

By attending to the formation of the successive expressions s*(Pj....P m ), it will readily 
be seen that the series 
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is a solution of the original equation, if 
M =l, 
M,=P.M , 
M S =P._ I M 1 +Q.M , 
M,=P._,M,+Q^ 1 M 1 +R.M M 



which last is a general relation determining the coefficient of any term from the 
coefficients of the n preceding terms ; or from the coefficients of all the preceding terms 
when the number of these preceding terms is less than n. This relation is the equa- 
tion of formation, and may be regarded as universal, bearing in mind that when p 
is less than n some of the terms of this relation vanish. 

The solution of the original equation is therefore reduced to that of a similar 
equation without second member : and, by what has preceded, this solution is 

M =f(P P ^ 

it being understood that M =l ; and this is the value of M^ before given. 
But the value of M p may be found otherwise, thus : 

Make P*- P +i=P^ Q*- P +2=Qpj &c. &c. 

Then the solution is evidently 

m,^(p;..p;), 

where the operation v has the same meaning as before, except that it is applied to 
the accented letters. 
Consequently 

where 

A p = 24 l q p+l + 2 2 r p+1 + . . -f- 2 /x ^. 1 z p+li 

ij p ::= Z 3 {q p+l A p _ 1 ) -f- Z 5 (v p+1 A p _ 2 ) +.. .-f- Zzx-iyZp+iAp-x+i) , 

and generally the terms are formed as stated in section 5, using p for x. 

9. I shall conclude this part of the subject with a few simple examples for the pur- 
pose of illustrating the processes here given. 

Ex. 1. Let the equation be 

That part of the solution which is independent of G*, is 

^=c{a*+(*~lX~ 2 &^ 



u. 



2 m 2 



U w : 
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and since in the present case P, and Q^ are constants, we have 

whence the particular value of u x is a series of terms of the form 

G w +aG^ l +(a 2 +b 2 )G^ 2 +(a 3 +2ab 2 )G^ 3 + 

The difference in character between the solution proposed in this paper, and that 
which would result from a perfect analysis of the general equation, may be exempli- 
fied by the present instance. 

The perfect solution of this equation is known to be 

where a and |3 are the roots of t 2 —at—h 2 z=.Q 3 and each 2 introduces a constant. 
Now if the constants be made equal, and the expression he written at length, we 
shall obtain the form derived above. 

The expansion of the two terms within the parenthesis gives 

G w+l +ttG x + u 2 G x „ 1 +..+a p+1 G p +.. + ca; v 
and if the difference of these be divided by a— /3, we get 

>y2 (22 mP+1 /2» + l „X Rx 

u *t a _(3 ^-lT-T a — 6 p "»•' • c ot — B 9 

which is the same in effect as the result which we have obtained. 
Ex. 2. Let the equation be 

U x = XU X __ j -J- Q U v „ 2 • 



JL lit/ 11 



«.=<*(*+ ] ){l +al* ^ +« 2 2*(^ 2"^) + 

r x I a 3 1 1 

=c|r(^+i)~ar^+^r(^^i)~^r(^^2)+..j; 

which may be put under the form of the definite integral 

u x =c j $r v ~°i;v x dv. 
If we add a second, member G w the equation for determining M p is 
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whence 



M,=^..(x- ?+ l)[l+a2? ( -_i^ 
=*...(*-p+l)|l+a( ir= i TT -j) 



1 

'•J 



+ « 2 fo(( 



1 1 \ 1/ 1 1 \\ ( 

V2V(*-iJ+2)(«-p+i)~(«-i)(«-2)y~*v«^Ti — *-2y/"r"" 

from which the series for u„ can be expressed. 
Ex. 3. Let the equation be 



Then 
u 



\ 



Zfte. 4. Let the equation be 



n 



U w A-lT'' ^#-2» 



# 



Then 



a? 



■h 



& 



\ 



■'~W+\)\} + ? 2 M*+ 1) + i* 2*(*(a?+ l)2- , a?(j?+ 1 )) + ••] 



+^(e (*+l)*(*~ l)(a?-2)(a?-3)(ar-4)+ | («+ l)x(x— \){x-2)(x-3)^+.. 



Ex. 5. Let the equation be 



Here 



„_ r XT' —1 TV" — *(*+*) (p + l)(p + «) 

» _ (g + l)g(g-l)(g-2) (p + 2)(p + l) (g+l)a (p + 4)(p + 3)(p + 2)(p + l) 
•^Vp — 2.4 2 ' 2 "I" 2.4 ' 



(a? + l)a?(a?-l)(a?-2) (a? -3) (a? - 4) (p + 2)(p + l) (x + l)x{x—l)(x — 2) 

__ 



*>p 2.4.6 



2 



N 



(m+2) 



■ (p + 4)(p + 3)(j> + 2)(p+l) (^ + l)^ _ (^ + 6)(^ + 5)(^ + 4)(^ + 3)(jp + 2)(^ + I) 
"^ 2.4 * 2 2.4.6 * 

(#+!)..(#- 2ro) _ Qp + 2)(p + l) (# + l)..(g-2ro + 2) (j> + 4)..(j? + l) (# + l)..(a?-2m + 4) 
2.4.. (2m + 2) 2 * 2.4. .2m + 2.4 * 2.4..(2m~2) 

, (p + 2m)..(jo + l) (a? + l)a?- r (p + 2m + 2)..(p + l) 
^'*— 2.4. .2m 2 ~ 2A..{ c 2m + 2) 
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The particular solution, therefore, is 

+ (3^ 2 ~5^-2)(P A ,..P,_JG,_ 5 +^ 
and the part of the general solution found by this method is 

U x — C{r x . . r ) i I -p 2 1 2.4 ' 2.4.6 ~T" , **J* 

Ex. 6. Let the equation be 

^xi®\*x^x-\\™ s 2^x*x-\Ux-<L\ % •\tt n \JL#..± x „ n+ i)Ux__ n 'ZZZvJ x ; 

then the equation for determining M^ is 

or making 

P -— P' 

* x—p + 1 ■*■ p) 
■M-p I &i±pWl p -i -J- G2*p*p~l Wl p __2 ""» • • r ^n\^p • • *p-n+ 1) N±p-n' =:Z ^ 5 

the solution of which is evidently the solution of 

M^+«iM^_ 1 +a 2 M^ 2 +..+a w M i? „ % ==0, 
multiplied by (P^..Pi). 

Let (^+fli^ l +o»^" a +..+^)" 1 =7^+f^ + ..., 

then M^=(P^.P' 1 ){c 1 A^ +w +c 2 Bi3^+..}; 

and taking the parts affected by each constant separately, it will be seen that the 
original equation reduces itself to a set of equations of the first order, 

v x - uV'jD^—cT^AGv 

v M -.pPj>._ 1 =p- 1 BG g9 



so that its complete solution can be adequately represented. 

Solution of Differential Equations in Series. 

10. I now proceed to point out a method by which the processes above indicated 
may be made to give solutions of certain general forms of linear differential equa- 
tions. 

In a paper on Linear Differential Equations presented by me to the Royal Society, 
and which the Society has done me the honour to publish in the Philosophical Trans- 
actions (Part I. for 1848, p. 31), I have enunciated, and so far as is material to the 
present purpose, demonstrated the following theorem : — 

That if, in a linear differential expression <p{x 9 D)^=X and its solution u=z^(x, D)X, 
the letter x be changed into the operative symbol D and D into —x 9 we shall thus 
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obtain another linear differential expression <p(D, — #)w=X, the solution of which will 
be u=-<p(D,-- #)X. 

In the application of this theorem, care must be taken that the first-mentioned 
solution is so written that the operations included under the function %|/ are not sup- 
pressed ; and it must also be borne in mind that the expressions obtained by this in- 
terchange of symbols will not in all cases be obviously interpretable. 

For applications of this singular analytical process I beg leave to refer to the me- 
moir above cited, where it is employed for the solution, in finite terms, of extensive 
classes of linear differential equations, and equations in finite differences. So far as 
the process is legitimate, it is to be observed that it is founded on reasoning of a 
purely analytical character. It does not in any manner whatever flow from the 
calculus of operations, or depend for its validity upon the soundness of the logical 
basis on which this calculus rests. 

Now it is a remarkable property of this mechanical interchange of symbols, that it 
instantaneously converts a linear equation in finite differences into a linear differen- 
tial equation; so that wherever the former is soluble, the latter is soluble also, pro- 
vided the result be intelligible, a condition always satisfied when the functions em- 
ployed are rational algebraical functions. 

As an instance worthy of notice, let us take the example last above given (Ex. 6.). 
Bearing in mind that u^ n =e~ nJ) u ay the proposed interchange gives the equation (writing 
<px for P*), 

tH-«iP(D)(rtO+^(DM^ 
whose solution, therefore, depends upon that of 

v-u<[)(D)(fv)=:ci n - l AG, 

a proposition established by Mr. Boole by the methods of the Calculus of Opera- 
tions. 

1 propose, therefore, now to employ this theorem of the interchange of symbols for 
the purpose of converting the forms of solution, above given, of equations in finite 
differences into the particular solutions of some general forms of differential equa- 
tions ; viz. those equations whose factors do not contain any irrational or transcen- 
dental functions of w, or contain them only in the form of series of ascending powers 

of x. 

11. Mr. Boole, in his General Method in Analysis, has shown that expressions of 

this character may be placed in the form 

/ (D)^+/ 1 (D)(^)+/ 2 (D)(a 2 ^) + ..=U, 
by changing the independent variable from x to its logarithm 6, and making use of 
the relation, ( D being -™ 



J){B-\)..{D-n+\)u-x n (J^ n u. 
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Making use of this relation, we immediately convert our equation, which we assume 
to be in the form 

d n u d n ~^u 

+ ..., + (a +b x-{- . . + k$x p ~ l +Z # P + . .)u= G, 
into 

« a (D..(D-»+l))«+ i w ((D-l).-(D-»))(i«tt)+ c n ((D-2)^(D-»-l))(i^tt)+... + W(D-/?)...(D-»- J p+l))( g Ptfa)+„ ^ 

+^-i((D-l).. (D-n+l))(«M+*»-i((D-2).. (D-» ))(s 2 ^)+....+^_ 1 ((D- J p).. (n-n-^))^*)*.. 

+« w _ a ((D-2).„(D-n+l))(.^t^ 

+ +.. 

or 

V+ ( D"— J (S * ) + \ 555=15 — ;/(■*") 

Now if in this equation we change D into 6 and 6 into — D 5 we obtain the equation 
in finite differences (which suppose to be of the nth order), 

a»w ff H — L — j -%-i + ^^j u e ^+....=(d(d~~l)..(d-n+l))~ 1 G d . 



>=s w *G«„...(2.) 



-«? 



or 



tttf+^w^ 1 +^«i^ 2 +^jM^8+.-..==(/o^) 'G^—H,^ (suppose); 
the solution of which, by section 8, is of the form 

where M =l, 

M,+f^M,+|M.=0, 



iV1 »+/ (fl««i+l) iV1 «-i"r/ (j« OT + 2) Mw ™ 2+,,, + /7F^ M ^ =0 ' 

Restoring the symbols, and thereby converting H*_ n into (/ (D))- 1 (g wfl G), which call 

H, we have 

^=M H+M 1 (^H) + M a (^H) + .. + M M (s^H) + ..^ 

where M , M x ... denote a series of operations having the following significations and 
relations ;— 

JL V JL a — •— A » 

M...i jjj~j 1VT n 
1 ~^" /* (D) 0~"~" * 



FINITE DIFFERENCES AND LINEAR DIFFERENTIAL EQUATIONS. 273 

M2 +/ (D-1) Ml +,/o(D) M ° = ' 



M \ M' D ~ m + 1 ) M | /a( D - OT + 2) / r (D-ro+r) -- _ 

mi "/o(D-m+l) llm -'" t "/ (D-m+2) iV, »-2~«~ •••+ f Q (j)„ m+r ) J-Vl m _ r _u ; 

or, by passing s s , s 29 , &c. outside the operations by the equation 

9(P)(f"H) =i"*p(D+»»)H, 
M=M H+g"M 1 H+s 2< 'M 2 H+..+2 m9 M m H+..., (3.) 

where the general law of relation is 

Now the expression H, which is the subject of all the operations, contains n arbitrary 
constants, since / (D) is of the nth order ; that is, provided a n is not zero. 
Let 0! j3 2 ...|3 n be the roots of f t=Q ; then the complete value of H is 

(/ (D))- l (tf*G) + Cl &+cJ*+ .. +cjf^ ; 

of which the first term will give us the particular solution of the original equation, 
and each of the other terms a complementary solution. Bearing in mind the equa- 
tion (p(D)(z p9 )=(pp.e p9 , we see at once that the first complementary solution is 

c/^(A +A/+ A/ + .. + A m e me + ..), 

or 

c 1 «£ /3l (A +A 1 <r+A 2l x 2 + # # ^ A m «r w + ..), 

where A =l, and the law of formation of the coefficients is 

and the remaining complementary solutions merely require the substitution of /3 3 ..*/3* 
successively for (3 15 with new constants. 

The reader will not fail to perceive the peculiarity of the series when a n is unity, 
(which value it can always have when it is not zero) ; for in that case the roots of 
f§t=Q are the natural numbers from to n— 1 inclusive, so that the series begin re- 
spectively with the terms 1, x 9 x 2 , ...a?*" 1 . 

If a n be zero, that is, if the factor of the highest differential coefficient of u do not 
contain an absolute term, then in order that the transformed equation may take the 
form (2.), it will be necessary to pass £ outside the operative functions in each term, 
and divide by s e . The initial function/ £ will then be of the form 

b n (t...(t—n+l))+a n „ x (t...(t — n+2)) 9 
and H will be 

C/o(D))-*( s <-*GV 

Similarly, if in addition to a n being zero, we have also b n and a n ^ respectively equal 
to zero, it will be necessary to pass e 26 outside the operative functions, and divide by 

MDCCCL. 2 N 
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£ 2d ; and so on if other terms of the factors are wanting. Thus we can in all cases 
obtain the transformed equation in the required shape ; but in all the eases where 
coefficients thus vanish, there is the important qualification that f t is no longer 
necessarily of the rath order ; so that H does not necessarily contain the proper 
number of arbitrary constants. The consequences of this consideration will be after- 
wards developed ; in the meantime we proceed to consider what modifications the 
series undergo where f t has two or more equal roots. 

12. If there should be two roots off t=0 equal to (Z 19 one series will be deficient ; 
and it will be supplied as follows. The expression H, or (/ (D))~ 1 (e^G) > contains in 
that case a term of the form k x fa hB ; and it is easily seen that 

p(D) (&*') = ty(D) (8**) + <p f (D) (**') . 

The wanting series is, therefore, 

k^Qog x( 1 + A 1< r+ A 2 # 2 + . .) + (Ai#+ A' 2 x 2 -\- ..) }, 



where 



Al= Ifl 



If there should be three roots of/ £=;0 equal tof3 u two series will be deficient, one 
of which will be supplied as last mentioned ; and the other by the introduction of the 
series 

& 2 <^ l {(log#) 2 (l+A 1 #+ 

where A'^-^i &c. ; for we have a term kjh^ ; and it is easily seen that 

9(D)(^¥^=^(D)s^+2^(D)s^+^(D)s^. 
And generally since, where there arep+1 equal roots /3 15 we have terms 

and since 

the deficient series will be supplied by the following, taken with a different constant 
for every value of p from unity upwards: 

&^j(log#f(l^ 

+^^(log<^-^ 

As a matter of convenience, when the equal roots are zero, a temporary nominal 
value should be given to them for the purpose of differentiating A n &c. 

13. Hitherto we have attended especially to the complementary solutions, or in 
other words, regarded G as zero. The operations, however, indicated in (3.) may be 
readily performed when G is a rational function of x ; which we will suppose to be 
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cleared of fractions, and so cleared to consist of a set of terms u p tV p . Then for this 
term C/o(D))- 1 (^G) =«,(/ (D))-V« + ^==a^(/ (»+ P ))- V« + ^ ; 

and we have, as before, 

u=cc p x (n ^\B +B l x+B 2 x 2 +..+B m x m +..), 

where B =:l, and the law of formation is 

fo(n+p+m)B m +f 1 (n+p+m)B m ^ 1 +f 2 (n+ 

and the whole solution is found by taking all the values of p. This will undergo a 
slight alteration, as the incipient term will be of the formf Q (D)(e (n ~ r)e G), where r factors 
of (2.) vanish by reason of some of the constant coefficients a n9 &c. being zero. 

In like manner it would be easy to represent the series if G contained log x and 
its powers ; but for most other forms it would be necessary to expand G in order 
to represent the series explicitly. The solution however is theoretically complete, 
since it consists solely in the performance of operations w T hich are known explicit 
functions of D. 

14. Before proceeding further with the main subject, I shall illustrate this process 
by a few examples. 

Ex. 1. Let the equation be 

(l+fc 2 <r+c 2 ^ 2 )^+K+V)^+« ^=G. 

Referring to (2.), we have 

/ D=D(D-l),/ 1 D=ft 2 (D-l)(D-2)+« 1 (D-l), 

/ 2 (D)=c 2 (D-2>(D-3)+6 1 (D-2)+a ;. 

roots of/^=0 are and 1. 

First complementary series^ 

c l (A +A l x+A 2 x 2 +..+A m x m -\-..), 
where 

A = 1, and m(m- 1) A m + (6 2 (m~ 1) (m^ 

is the law of formation. 

Second complementary series, 

c 2 ^(A +A 1 o?+A 2 ^ 2 +..+A^ m +..), 
where 

A =l,and(m+l)mA m +(& 2 (m)(m— l)+a l m)A m ^ 1 +(c 2 (m^ l)(m—2)+b l (m—l)+a )A m ^^O 

is the law of formation. 
Particular solution, 

where 

H =(£r G >^o(D)=^ 

is the law of the formation of the operative functions. 

2 n 2 
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The three laws of formation, which are here written down at length, are in sub- 
stance the same ; for from the last the others are made by merely writing the roots 
successively in lieu of D in the factors. 

Ex. 2. Let the equation be 

ujU du 

(rf+qa?) j-2 + (x+px 2 +5qx*) s + (px—n 2 +(4q+r)x 2 )u=G. 

Referring to (2.), we have in the first instance, 

((D-2) 2 -ra 2 )(s 2 ^)+jp(D^ 
whence 

(D 2 -n 2 )u+pI)(z 9 u) + (qB 2 +r) (e e u) =s~ 2 *G ; 

roots off Q t=0 are n and —n. 

First complementary series, 

c 1 #%A +A 1 #+A 2 r 5 + .. + A m x m -+- . •)> 
where 

A =l, and ((w+m) 2 ~w 2 )A m +^(^+m)A^ 1 +(9(n+m) 2 +r)A OT „ 2 ==^ 

is the law of formation. 

Second complementary series, 

where the law is, as before, changing the sign of n. 
Particular solution, 

M=^ (D)H+^ 1 (D)H+..+g^ B (D)H+.., 
where 

H=(D 2 ~<r(s~ 2 *G), ^ (D) = 1, 
and 

is the law of the formation of the operative Amotions. 
If p=0, 9=0, and r=l, so that the equation becomes 

d u du 

the law is 

A — * A. m —2 

m ~" M ' m(2n + m) 

whence 

+ c«*-"{l-i(l-*)- 1 ^^ 
Ex. 3. Let the equation be 
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Referring to (2.), we have 

D 3 u+q? $ u=G; 

/ (D)=D 3 ; three roots equal to zero; 

/.(D)=y; 

m 3 A m +§ , A ro _ K =0, law of formation. 
Hence the complementary series are, 

"" Cl [ 1 "^ itt rf(ai)»* ~" n 3 (2n)« (3w) 3 * I*"" 

+d V«« 3 V»^W« 8 (2w) s + V w + 2n ^ 3m/ n 3 (2m) 8 (3m) 3 * "•""V' 

-f C 3 |(l0g«) 2 ^l - tfX n + ^J2nf^ n — ^(2^)» (3^8 **"+•• ) 



/ 



+ 6 log X (^ n3 #" — ^ + 2n j w3 ^ 8 * 2W + ^ + 2n 4- 3n J ft8 (2m)3 ^p 3**+ . . j 

■ r/ 2< ? r"4-f 2 4- 2 _|_ 3 V g r ». 

"f" ^ - n 2 n 3,r "^Vw 2 ( 2w ) 2 ri{2n)Jn 3 (2n) aX 

/I. _? l 2 . 3 . 3 . 3 Y£._? 9_^, \1 

~ V« 2 T" (2m) 2 "T" (3m) 2 "T n(2n) "*"n(3») " l "(2w)(3«)y« 3 (2m) 8 (3m) 3 * "^ "yj ; 

and the particular solution is 

M=H-yx re (D+M)- 3 H+ ? V M (D+«)- 8 (D+2«)- 3 H-..., 
H being B~ 3 (f B G)- 

15. The completeness of the preceding forms of solution depends, as above inti- 
mated, upon the circumstance that the function f (D) is of an order not lower than 
the order of the original equation. It may however be of a lower order, as would 
take place in the first of the examples above given, if a 2 instead of being 1 were zero 
and b 2 were also zero. 

Let the equation then be {Ex. 4.), 

d u u/ij 

c 2* 2 &*+ + V)^+«o M =G. 

Referring to (2.), we have, in the first instance, 

(D-l)(^)+(c 2 (D-2)(D-3)4-& 1 (D-2)4-a )( 6 ^)=G, 



or 



Dw+(c 2 (D-.l)(D~2)+5 1 (D~-l)+«o)(^)=^G 
so that/ (D) is of the first order, and the root is zero. 
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The only complementary function to be obtained is therefore 

f l l 

w=c 1 ll — a x^a (b l +a )^ 2 ^^a (b l +a )(l.2c 2 +2b l ^a )o^ 

+2lA a ^ +a ^ l ' 2c ^ 2bl+a ^ 

a series which is divergent, but which, as will also be seen afterwards, is finite when 
the constant coefficients are connected by the formula (p — • l)pc 2 -\-pb l -\-a i i=0, p being 
any positive integer. 

Further, it may happen that the series obtained by the process, even when they 
afford a complete solution in respect of the number of constants, are divergent for all 
or for some values of at. This may evidently be the case in the solution of the second 
of the examples above given ; for the law of the coefficients, as we advance in the 
series, approximates to A m = — #A TO _ 2 . 

In these cases, other solutions in series may be obtained by resolving the equation 
in finite differences in a series of terms of the form H e+P instead of H d „ p . In order 
to effect this, all that is necessary is to write 0+r for 0, (r being the order of the 
equation in finite differences), and to divide by the factor of the last term instead of 
the factor of the first term ; or in other words, we must pass ? 9 outside the functions 
in (2.), and multiply by e~ rd ; so that this equation now assumes the form 

and the equation in finite differences is 

We have now to inquire for the roots off r t=0 ; the incipient term is (jC(D))-' 1 (e (n ~ r) *G), 
which call H! ; and the particular solution will be found to be 

where 

4o(D)-l,and/^D^m)^D)4/^ 1 (D^m)^ 1 (D)+...=0 

is the law of the formation of the operative functions. 

The substitution of the roots of f r t successively for D gives the law of formation of 
the complementary series. 

Taking the last example, the transformed equation now becomes 

(c 2 D(D^l)+& 1 D+a o >+(D+l)(g- w)=g- 20 G. 

Let the roots of f r t or c 2 t(t— l)-\-b 1 t+a Q =Q, be /3 X and j3 9 . 
The first complementary series is 

c 1 ^(A +A 1 x~ 1 -\-,,+A m x- m +..), 
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where A =l, and/ r (|8 1 — m)A m +((3i— '^+l)A m _ 1 =0 is the law of formation ; the series 
then are 

r /*&_ & jfii-*A- & l Pi~\... *&-*_ A- £lZ2 & ~ 2 yft-3 > 

C T /r(fc-lf ^(fr-l)/^-^ /r08i-l)/r(^,-2)y;03 1 -3) a? +•* 

+ c 2 { similar function of |3 2 '.........}, 

which is, (since/^jS!— w)=»(w+l)c 2 — SwjSj—wi!), 

^ J ^ & rA-ig. & & — 1 r g,-2 

Cl [* 2c 2 -2/3 1 -6 1 tX "T"2c 2 -2^ 1 -* 1 2.3^-4/3^26/ 

ft ft-1 ft -2 ft ] 

2c 2 -2ft-6 1 2.3c 8 -4/3 1 -24 1 3.4^-6^-36/ ' +**/ 

+ c 2 { similar function of |3 2 ...........}. 

When (3 X or j8 2 is a positive integer, one of these series is terminable ; and if both 
are positive integers, the series derived from the smaller root is terminable^ and the 
other gives no result, the coefficients becoming infinite. The first of the series will 
then be found to give the same result as that produced from the divergent series 
(which is evidently terminable in form in the case indicated), except that it begins at 
the other end. In this case the other complementary solution can be found in finite 
terms by reducing the order of the equation. 

The particular solution is 

u— n x x ^D^^rii-t-^ y>(D-i)/ r (D--2) "i"""— > 
H, being (/.(D))- 1 (g- 2 ^G). 

Let us now return to example 2, the solution of which, as above found, is in some 
cases divergent. 

The transformed equation now becomes 

( ? (D+2)^r)w+j9(D+2)(r%) + ((D+2) 2 ^r 2 )(r%)=s^G. 
Roots of/i—0 are —2+ \f (~\ which call ^ and j3 2 . 
First complementary series, 

where 

A =l, and (q(^ 1 -m+2) 2 +r)A m +p((i 1 — m +2)A m _ l + ((p i — m+2) 2 ~r 2 )A^ 2 =0, 

is the law of formation ; and the second complementary series is the same, using the 
other root with a different constant. 

In like manner the particular solution is easily represented. 

16. We are now in a position to discuss the character of the various series obtained 
by this process with reference to their convergency or divergency, a subject of the 
highest importance to the value of the process. The investigations which follow, 
will, it is apprehended, be found to afford a complete test of the nature of the series 
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In general, it will of course be understood that these researches, as to convergency 
and divergency, relate only to the complementary series, or in other words, to equa- 
tions deprived of their second member; nevertheless they will to a great extent, if 
not throughout, apply to cases in which the second member is in the form of integer 
powers of x or of log x. 

1st. When in the set of functions 

/.(D), /(D), /(D), .../_,(!>), /(D), 

the function f (D) is of a higher dimension with regard to D than any other of the 
set, or is what we shall here call the dominant function, the solution of the equation 
can always be found in a convergent series of ascending powers of x ; and if in such 
a case we solve the equation in a series of descending powers of x 9 which we can do 
if we please, that series is certainly divergent. 

This is immediately apparent from the consideration of the law of the coefficients 

which, as m increases without limit, approaches to the form 

A — _^A — *^A — — ^A nvm\ ——n A — h A — 

^m — /"m f m m ~~ 2 * * f m , m ~ 1 " f "i^m""" Ll n—l jrx m—l u n—l- cx m—2 ••••5 

assuming that all the functions are only one degree lower than A m9 which is the least 
favourable case for convergency. Therefore, if the largest of the terms of the right- 
hand side of this equation be 4_ 1 A m _ i0 , we have 

A < rln - l \ 

and we can therefore arrive at a point in the series at which the ratio of the coeffi- 
cient of x m to that of x m ~ p diminishes without limit. 

It will also be observed that the series introduced by two or more equal roots are 
of the same character as the original series from which they are derived ; for, A m being 

dA. dA 
of the form <p(^i-\-m) 9 we have *^ =r ~^f » ail d when A m „ p is of a higher order with re- 

-p^ 2 is also of a higher order than -j 1 

dm ° am 

ferential coefficients. 

We have now merely to inquire what must be the form of the original equation 
t,hat/ (D) may be the dominant function. Referring to (2.), we see that it is neces- 
sary that the factor of the highest differential coefficient of u should contain one term 
only. If this factor be 1 or x 9 no restriction need be imposed on the succeeding 
factors. If it be x 2 , the factor of the next lower differential coefficient must not 

d n u d n ~*u 

contain an absolute term ; and generally, if x p be the factor of j-^ the factor of ■ , n-l 

must begin with a term not lower than x p ~\ that of ^=2 w *th a term not lower than 
x p ~ 2 9 and so on. 



ference to m than A m , 7 - is also of a higher order than -j-^, and so for the other dif- 
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In short, all equations where w does not enter into the factor of ^> and all equa- 
tions of the form 

where p is a positive integer, are soluble in convergent series of ascending powers of w. 
The third example above given is an instance of this form. 
2ndly. When in the set of functions 

/ r (D), ./UD), .../(D), /«(D), 

the function / r (D) is the dominant function, the solution of the equation can always 
be found in a convergent series of descending powers of x ; and if in such a case we 
solve the equation in a series of ascending powers of *r, that series is certainly diver- 
gent. 

This is apparent, as before, from the consideration of the law of the coefficients, 

which, as m increases without limit, approaches to 

in the least favourable case for convergency. 

On proceeding to inquire what must be the form of the original equation, we see 
again that it is necessary that the factor of the highest differential coefficient of u 
should contain one term only. If this be x p 9 then the other restrictions are, that the 
factor of the next differential coefficient must stop at x p ~\ that of the next at x p - 2 , 
and so on. 

In short, for all integer values of p the equation 

is soluble in a convergent series of descending powers of x. The 4th example above 
given is an instance of this form. 
3rdly. When in the set of functions 

/.(D), /(D), .../..(D), /(D), 

the functions / (D) and/,(D) are of the same dimensions, and are both dominant 
over all the other functions, the solution of the equation can be found in a series of 
ascending powers of x 9 which for some values of x is convergent, and for other values 
of x is divergent ; and the solution can also be found in a series of descending powers 
of x which is divergent for all values of x for which the other series is convergent, 
and convergent for all values of x for which the other series is divergent. 

For in the ascending series the law of the coefficients approaches the form, (the 

MDCCCL. 2 O 
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coefficient of the highest power of m mfjn being 1, and that in f r m being / n? ) 

A — — / A 

J - x m v n £X m~rJ 

and in the descending series the law approaches to 

/ a - — I- A 

The former series is therefore convergent for all values of x numerically less than 

(4) ~ 7 j an d divergent for all values of x numerically greater than this limit ; and the 

latter series is divergent for all values of x numerically greater than this quantity, 
and divergent for all values of x numerically less. 

The equations to which this rule is applicable are of the forms, (p' being less 
than p,) 

(^+^)^+(&-,aH'- , +..+^ 1 ^ 
and 

The second example above given is an instance of the first of these forms. 
4thly. When in the set of functions 

/.(D), /,(D),.../_,(D), /(D), 

one or more of the intermediate functions is or are of the same order as the extreme 
functions / (D) and jC(D), or as the highest of these two when they differ in dimen- 
sions, the series obtained by the above processes will be divergent for some values of 
x, and we have not as yet any method of deriving a convergent series corresponding 
to these values ; and if one or more of the intermediate functions be of a higher 
dimension than the extreme functions, the series obtained by the above processes will 
certainly be divergent. 

These remaining cases therefore sever into two species ; first, where some of the 
intermediate functions are of the same order as the highest of the extreme functions; 
secondly, where one or more of the intermediate functions are dominant. 

The first of these species includes equations of the two following forms : — 

d n u d n ~ l u 

in which there can be no function higher than/ (D) ; and 

d n u d n ~ l u 

in which there can be no function higher than f r (D). 

In these cases it will easily be seen that the law of the coefficients of the ascending 
series, as m increases without limit, approximates to 
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and the series therefore approaches without limit to a recurring series, in which the 
constants of relation are 

___ nr* mmmmm rp& mmmm ^ __ /yj/' • 

iAj , <Aj , • • • *Aj 9 

and the denominator of the rational fraction, to which the residue of the series ap- 
proaches, is 

In the descending series, the law of the coefficients approximates to 

and the series approaches without limit to a recurring series, in which the constants 
of relation are 

k n 1 u n 1 u n 1 t 

t/fl X Iff, vC" v<yi vC" 



and the denominator of the rational fraction, to which the residue of the series ap- 
proaches, is 

When/ (D) is higher than/ r (D), the ascending series alone can be used ; when/ r (D) 
is higher than f (D), the descending series alone can be used ; and when f (D) and 
f r (D) are of the same dimension, either may be used ; and the approximations above 
referred to render it probable that these series, notwithstanding that they may be 
divergent, are the developments of continuous algebraical expressions. 

The second of the species above referred to includes all equations which are ex- 
cluded from the preceding forms ; that is, all forms which transgress both the restric- 
tions to which the equation in the third case is subjected. 

Of these forms, the solutions, whether obtained in ascending or in descending series, 
are always divergent ; and the divergency appears to be of an extreme and unmanage- 
able character. In this case we have an intermediate dominant function ; and the 
convergent solutions might, from considerations of analogy, be presumed to be series 
infinite in both directions, the roots of the dominant function determining the inci- 
pient terms. 

The treatment of these forms requires the solution of the equation in finite dif- 
ferences 

..+RX +3 +QX +2 +PX +1 +t^ 

not starting from either of the two extreme terms, as is done above, but from the 
term u x9 so as to get a result in the form 

»,= ..+M. J G l+ i+M. 1 6 #+ i+M 6 l +M 1 6 l . 1 +M,6 M +.. 

It would probably not be difficult to show that such a solution exists ; but I have 
not found one in a form available for the purpose to which it is desired to be applied. 

2 o2 
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17. Throughout the preceding investigations, the series obtained by the processes 
here displayed undergo a modification of form in the event of the expression f t or 
f r t, as the case may be, having one or more sets of imaginary roots. 

Let there be a couple of such roots of the form a+|8\/— 1. In the ascending 
series these roots give 

c 1 x« + ^' r[ (A + A,x+ A#' 2 + . . +A m x m + . .) 

c 2 x-^~(B +B 1 x+B 2 x*+ . . +B m x m + . .), 



where 



and 



/o(«+^+f3 N /-l)A m +/ 1 (a+m+f3 N /- 1 ) A m~i+...=0 



/•(«+w-^-l)B lli +/ 1 (a+i»-^ x /-l)B III _ 1 +...=0 

are the respective laws of formation. 

It is apparent, therefore, that if A m be of the form ^(a+|3^/— 1), B m is of the form 

Making c 2 =c u and remembering that 

x ^~i +x -^~i =2 cos (|8 toga?) 
a*^ 1 — ar^~=2^^sin (filogx), 
the sum of the two series gives the double series, 

• 2cv*"(cos (]31og^){A +B +(A 1 +B 1 >+.. + (A m +B m K+..} 

+ N /~lsin(/31og^){A ^B +(A 1 -BO^+..+(A m -B m )^+..}) ; 

which is necessarily real, since A m +B m is purely real, and A m — B m is purely imaginary. 
Making now c l = — c 2 , the sum of the two series gives another double series, (making 

c a = — A\/-l,) 

2^ a (x/^cos (j3 log *0 { A ^ 

+ sin( i 31og^{A +B +(A 1 +B>+.. + (A w +BJ^+..}), 

which is likewise real. 

The descending series may be treated in a similar manner. 

18. Most of the examples to which the preceding processes are applied have been 
taken from the paper in the Philosophical Transactions for 1844, in which Mr. Boole 
developed his new General Method in Analysis, with which the subject matter of the 
present paper is closely connected, though the methods exhibited are distinct ; unless 
indeed it should prove, that the interchange of the symbol of operation and the inde- 
pendent variable, and the general relation exhibited by Mr. Boole's fundamental 
theorem of development connecting any system of linear differential equations with 
a corresponding system of equations in finite differences, are merely different repre- 
sentations of a part of some more general method or process. 

The principal difference in results, so far as concerns the solution in series of linear 
differential equations, appears to be, that in this paper the law of relation of the 
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coefficients of each series is distinct, and substantially the same in form for all; and 
that it is not necessary to have recourse to the method of parameters in the cases of 
equal roots, or in the case of there being a term G on the right-hand side of the 
equation. In the case of imaginary roots, the laws of the series are in the first instance 
distinct from each other, and afterwards combined in couples. 

19. The investigations contained in the latter part of this paper reduce the problem 
of the integration in finite terms of the general linear differential equation with 
rational coefficients, to the finding of an algebraical expression representing the deve- 
lopment 

A I A np \ A /y>2 I !._ A /y>Wi t 

where A =l, and the law of relation is 

/ ^A m 4/ 1 mA m _ 1 +/ 2 mA m „ 2 + .. +/ r mA w ^=0 ; 

the functions fofif 2 --fr being known specific functions. The series to be summed 
closely resembles a recurring series ; it differs from it in this particular, that the law 
of relation, instead of being constant, has a uniform and simple variation as it pro- 
gresses along the series. If the rational coefficients should be themselves infinite 
series, the process still applies, the only difference being that each term would be 
formed from all the preceding terms, instead of being formed from the r immediately 
preceding terms, or from all the preceding terms when the number of them is less 
than r. 

In those cases in which the equation is soluble in finite terms by known methods, 
we are enabled to assign the algebraical expression for the series ; a result which may 
be used for the discovery of generating functions. 

Thus, taking the general equation of the first order, 

du 

we see that 

where A =l, and the general law of the series is 

a 1 mA m +(b l (m—l)+a )A m ^ x + (c 1 (m— 2) + b )A m „ 2 + ...(l 1 (m-n)+k )A m „ n =0. 
If we take the general equation of the second order, 

a a u/ii 

(a 2 +b 2 x+ .. +l 2 x n ) j^ +(a l +b l x+ . . +* 1 tf l ~ 1 ) ^ + (a +b x+ ..+h x n - 2 )u=0, 

the law of the series will be 

ajn(m— l)A m +(b 2 (m--l)(m-- 2)+a 1 (m— l))A w ^ 1 + (c 2 (m— 2)(?n—3) + b l (m--2)+a )A m ^+., 

+ (L 2 (m— n)(m— ^— l)+^(m— w)+^ )A m ^ w =0 ; 

and generally, in equations of the nth order, the law of relation involves the nth 
power of m, the number of the term sought for. Thus the determination of soluble 
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cases of linear equations resolves itself into the inquiry, in what cases can a series 
whose scale is of the nth order be resolved into a number of series whose scale is of 
the first order. 

20. I have not thought it necessary here to extend these solutions in series to 
linear partial differential equations. The process by which the extension can be 
made is well known, and has no peculiar relation to the methods here developed. 



